Nonlinear Spectral Singularities for Confined Nonlinearities 
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We introduce a notion of spectral singularity that applies for a general class of nonlinear 
Schrodinger operators involving a confined nonlinearity. The presence of the nonlinearity does not 
break the parity-reflection symmetry of spectral singularities but makes them amplitude-dependent. 
Nonlinear spectral singularities are, therefore, associated with a resonance effect that produces am- 
plified waves with a specific amplitude-wavelength profile. We explore the consequences of this 
phenomenon for a complex 5-function potential that is subject to a general confined nonlinearity. 
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Introduction. — A spectral singularity of a complex scat- 
tering potential is a mathematical concept introduced 
and studied by mathematicians for more than half a cen- 
tury [l|, Q ■ This concept that applies only for non-real 
scattering potentials, was recently shown to have an in- 
triguing physical interpretation |3| : It corresponds to the 
energy of a scattering state whose reflection and trans- 
mission coefficients diverge. This observation has moti- 
vated identifying spectral singularities with certain zero- 
width resonances and led to the study of their physical 
implications In optics, a spectral singularity gives 

rise to lasing at the threshold gain Q. Its time-reversal 
corresponds to a coherent perfect absorption (CPA) of 
light that is also called anti-lasing. The observation 
of the latter reported in Q may be considered as an ex- 
perimental evidence for the physical relevance of spectral 
singularities. 

Once one creates a spectral singularity in an optically 
active material [H, H| , it begins amplifying the background 
noise and functions as a laser. One may argue that a real- 
istic treatment of this phenomenon should also take into 
account the nonlinearities associated with the production 
of high-intensity radiation in the gain region. This pro- 
vides our main physical motivation for the study of the 
meaning and behavior of spectral singularities for non- 
linear operators, a task that to the best of our knowl- 
edge has not been considered previously @. The search 
for an appropriate nonlinear extension of the notion of 
a spectral singularity is particularly important not only 
because it has been an open mathematical problem for 
several decades, but also because it can pave the way for 
the discovery of the nonlinear analogs of the interesting 
physical phenomena such as threshold lasing [f| , antilas- 
ing Q, and CPA-lasers in areas other than optics. 

Motivated by the physical meaning of a spectral singu- 
larity of a linear Schrodinger operator 0], in this letter 
we propose a definition for a spectral singularity of the 
nonlinear Schrodinger operators 'H 1 of the form 

H 7 ip(x) := -tp"(x) + v(x)tp(x) +-fx(x)f(\ip(x)\,x)ip(x), 

(1) 

where v is a rapidly decaying complex scattering poten- 
tial, 7 is a nonzero real coupling constant, x( x ) := 1 f° r 
x 6 [0, 1] and x( x ) '■= for x £ [0, 1], and / is a real- 
valued function. The presence of x m © shows that 
involves a confined nonlinearity. A concrete exam- 



ple is the confined Kerr nonlinearity, with f{\ip{x)\, x) := 
\ip(x)\ 2 , that appears in the study of Bose-Einstein con- 
densates (lOj and has well-known applications in optics. 

The problem of introducing spectral singularities for 
nonlinear operators is plagued with severe mathematical 
difficulties associated with the proposing an appropriate 
definition for the spectrum and a suitable scattering the- 
ory for these operators. The simple idea of considering 
confined nonlinearities, that is mainly motivated by phys- 
ical considerations [11], allows for circumventing these 
difficulties. As we show below this idea plays a central 
role in our ability for defining a useful notion of a non- 
linear spectral singularity (NSS). 

The time-independent nonlinear Schrodinger equation 
corresponding to ([T]) is given by 



H 1 ijj{x) — k 2 tp(x), 



(2) 



where k is a complex number. It is easy to see that 
outside the interval [0,1], © coincides with the linear 
time-independent Schrodinger equation; 



ip"(x) + v(x)4>(x) = k 2 %l)(x). 



(3) 



Jost Solutions and Nonlinear Spectral Singularities. — 
Consider the linear operator H := — d 2 + v(x), whose 
continuous spectrum is [0,oo). The eigenvalue equa- 
tion for H, i.e., ([3]), admits the so-called Jost solu- 
tions ipk± that fulfil the asymptotic boundary conditions: 



i>k±{x) -> N±e 
plex numbers N± 



as x — > ±oo, for some nonzero com- 
In particular, 



lim [ip' k ±(x) T ikipk±(x)] 



o. 



(4) 



The Jost solutions, ipk+ and ip^- , are the scattering so- 
lutions of ([3]) corresponding to incident waves from the 
left and right, respectively. 

Definition 1: A positive real number k 2 is called a spec- 
tral singularity of H if tpk± are linearly-dependent 0], 
i.e., ip k - oc ip k +- 

The following is a simple consequence of this definition. 

Theorem 1: A positive real number k 2 is a spectral sin- 
gularity of H if and only if there is a solution ipk of J3J) 
such that lima; _> -too ipk{x) exist as nonzero complex 



2 



numbers, and ip^ satisfies 



lim tp' k (x) T ikipk{x) = 0. 

:— >±oc L 



(5) 



We introduce a notion of spectral singularity for the non- 
linear operator © by promoting Theorem 1 to a defini- 
tion. 

Definition 2: A positive real number k 2 is said to be a 
spectral singularity of i? 7 , if there is a solution ipk of © 
such that \im x - > .± 00 tl>k(x) exist as nonzero complex 
numbers, and ipk satisfies ([5]). 

In what follows we use the term NSS for a spectral sin- 
gularity of the nonlinear operator (JXJ) . 

In order to ensure that the physical interpretation 
of spectral singularities is left intact, we demand the 
existence of the Jost solutions of the nonlinear equa- 
tion © . We identify them with the solutions that satisfy 
the asymptotic boundary conditions © , denote them by 

t/j^2 , and keep using ipk± for the Jost solutions of the lin- 
ear equation ([3]). Moreover, because the nonlinearity is 

confined to [0, 1] and are continuously differentiablc, 



^}{x) =Vfc-(z) forx<0, 
^k+( x ) = 4>k+{x) for x > 1. 



In particular, 



^(o) = ^_(o), vi 7 _ ) '(o) = ^_(o) ) 

^+(1), ^'(1) = ^ + (1). 



(6) 



(7) 

(8) 



We can view Eqs. ([7]) and ([SJ as initial conditions for 
the differential equation ©. Solving the initial- value 

problem defined by ([2]) and © for x > gives ip^} on 
[0, oo). Similarly, solving the initial- value problem de- 
fined by © and © for x < 0, we find ^1+ on (-co, 0]. 

These together with © specify throughout R. Note, 
however, that this procedure works provided that the 
above initial-value problems have global solutions. In- 
deed, because the nonlinearity is confined to [0, 1], it suf- 
fices to make sure that they have solutions on [0, 1]. The 
spectral singularities of H 7 are given by the values of k 2 

(V) 

for which at least one of the Jost solutions ipl± satisfy 
the conditions listed in Definition 2. 

Potentials Vanishing Outside [0, 1]. — InRefs. [l2l.[l3j. we 
consider linear spectral singularities (LSSs) of potentials 
that vanish outside [0,1]. This allows for making more 
definitive statements about the behavior of these spectral 
singularities. The same holds for the NSSs. 

Suppose that v(x) = for x £ [0,1]. Then, 

Eqs. ((H]) and ©-(0) respectively take the form ^}(x) — 

N_e~ lkx for x < 0, i/)j?2(x) = N + e lkx for x > 1, and 



tf2 



(0) 
(1) 



(0) 
(1) 



-ikN-, 



ikN. 



+ ■ 



(9) 
(10) 



where N. 



+ 



N + e lk . 



In order to determine ip k 2(x) for 
x > 0, we solve the initial- value problem defined by © 
and © on [0,1]. The solution, that we denote by 
gives the value of tp k 2(x) for x G [0, 1], provided that it 
exists. For x > 1, the general solution of © is a linear 
combination of plane waves. This observation together 



with the requirement that tp k 
tiable at x = 1 give 



(7) 



is continuously differen- 

for x < 0, 
for < x < 1, 



2ifc iU1 X 1 - 



(11) 



Here is the solution of © on [0, 1] that satisfies 

e fe (0) = iV_ ! &(0) = -ifcJV_, (12) 
and F±(k) :— ££(1) ± Similarly, we obtain 



■G + ( fc )-e-^G_( fc ) f X<Q 



Cfc(^) 

jy- e ifc(a;-l) 



for X e [0,1], (13) 
for x > 1, 



where G±(k) := Cjt(0) ± ifeCfc(O) and is the solution of 
© on [0, 1] that fulfils the initial conditions: Cfe(l) = 
and Cfc(l) = ifciV + . 

Because 1/4+ an d V'i- respectively correspond to the 
scattering states with incident wave from the left and 
right, we can use (fTTj) and (|TU| to determine the left and 
right reflection and transmission amplitudes, R and T. 
This gives 



R 
R r 



G + (fc)' 
- 2lfc F+(fc) 



e 



V 



2it 



G + (k) ' 
2ike~ ik N. 
F-(k) 



(14) 



where the superscripts / and r stand for left and right, 
respectively. For 7 = 0, G± (resp. F±) are proportional 
to N + (resp. iV_), and the latter drops out of (|14l) . 

Having obtained the explicit form of ip k 2, we can i m ~ 
pose the condition that they yield a NSS. Demanding 
that V'fcj- satisfies ©, we find G+(k) = 0. Similarly, im- 
posing © on ipy} gives F_(k) = 0. Therefore, in view 
of (Q3J , NSSs can again be interpreted as the energies of 
certain zero- width resonances [3j . 

The condition that cither ipjrf or ip k 2 gives rise to a 
NSS is equivalent to demanding that £ and £ are solutions 
of the boundary-value problem defined on [0, 1] by © 
and the outgoing boundary conditions: 

ii(o) - ifcVfe(o) = 0, 4i(\) + ifcVfc(i) = 0. (15) 

It is easy to see that for potentials vanishing outside [0, 1] 
Eqs. (|15p arc equivalent to ©. In particular, they are in- 
variant under the parity {V) transformation: x — > 1 — x. 
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This is a manifestation of the fact that, similarly to the 
LSSs [l4[, the ^-symmetry of the boundary conditions 
(O or (fT5|) leads to an intrinsic "P-symmetry of NSSs. 
This means that once the parameters of the system are 
tuned to realize a spectral singularity, it will amplify the 
background noise and begin emitting radiation of the 
same wavelength from both ends. 

The main difference between LSSs and NSSs is that the 
resonance effect corresponding to the latter is intensity- 
dependent. The system amplifies an incident plane wave 
of negligible amplitude and produces outgoing waves of 
the same wavenumber k and a particular fc-dependent 
sizable amplitude. This intensity-dependent resonance 
effect may, for example, be used to devise a measurement 
scheme that determines the wavelength of an incident 
wave using the information about the intensity of the 
transmitted wave. 

Complex 8 -function potential. — Consider the potential 



v(x) = $5(x- a), 



(16) 



where 3 is a complex coupling constant and a £ (0,1). 
This potential supports a single LSS provided that 3 is 
purely imaginary [la ] . For this choice of v and the func- 
tion / given by f(hHx)\, x) = |-^(o;)| 2 , Eq. ([2]) admits 
analytic solutions [Tg|. For real values of 3, this model 
has a ppl ications in the study of Bose-Einstein conden- 
sates [13. See also (l8j . 

The study of the NSSs for the potential (|16|) requires 
solving ([T]), that for < x < a and a < x < 1 takes the 
form: ip" + k 2 ip = ^ f (\ip(x)\, x)ip(x) ■ This is equivalent 
to 

^(x) = Mx) + j[ 9(?,v)fMv)\,vMv)*u, ( 17 ) 

Jxq 

where ipo(x) is the general solution of ip" + k 2 ijj = 0, 
&(x, y) := sin[fc(.T — y)]/k is the Green's function for the 
latter equation, and xq £ (0, 1) is arbitrary. Using this 
equation to express the tp(y) appearing on its right-hand 
side in terms of tpo and and repeating this procedure, 
we can obtain a perturbative expansion for i/j(x) where 
7 serves as the perturbation parameter. In the following 
we consider a homogenous nonlincarity where / does not 
explicitly depend on x, i.e., / = f(\ip(x)\), and perform 
a first-order perturbative treatment of the NSSs of (fTBf . 

First, we recall that substituting (fT5|) in ^ is equiv- 
alent to demanding that ip is continuous at x = a and 
that ip'(a + ) = tp'(a~) + $ip(a), where ip'{ a+ ^) stands 
for the right /left derivative of ip at x = a. Next, we use 
(|17|) to obtain a perturbative expression for the solution 
£fc of ([2]) in the interval [0,a) and use the continuity of 
ip at x = and the above matching condition for ip' to 
extend it to [a, 1]. Finally, we demand that the result 
also satisfies the second equation in (fT5|) . This gives a 
pair of equations that we can solve to express 3 and £fc(l) 
in terms of k,a, and iV_ := £fc(0). The result is 
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2ik 1 



lf-A 
4fc 2 



<?(7 2 ), 



(18) 



and & (1) = e~ 2lka N^ (l + + 0(7 2 ), where /_ 

f(\N-\), A := e ^k(l~a) + e 2ika _ 2^ and B ,_ e 2ife(l-a) 



-2ik(2a— 1). Eq. (fT8|) is the condition under which 

ipk- yields a NSS. For 7 = 0, it reduces to 3 = 2ik which 
determines the corresponding LSS [TJI, [l5| . Similarly we 

find that ipk+ gi ves r i se to a NSS provided that we enforce 
(fI5]l after replacing /_ with f + :=f(\N + \). 

Let r and s denote the real and imaginary parts of 3, 
so that 3 = r + is. Noting that /_ is real, we can solve 
(fT5|) for 7/_ and s in terms of a, k and r. This gives 



if- 



kr 



2k- 



sin£;cos[£;(l — 2a)] ' 
cosfccos[/c(l — 2a)] — 
sinA;cos[fc(l — 2a)] 



(19) 
(20) 



where we use ~ to mean that we ignore C(7 2 ). Because 
cosine is an even function, these equations are invariant 
under the 7- > -transformation: a — > 1 — a. Therefore, we 
can confine our attention to the case that a < |. 

Eqs. (fT9"f and (f2T)]) provide a reliable description of the 
NSSs of (HBJ provided that the right-hand side of (19} 
is much smaller than k 2 . This implies that < \r\ <C 
I sin/ccos[/c(l — 2a)]|fc. In particular, 



< \r\ < k, 



and for all integers m, 



7rm for all a, 

k + \ 7r(m + l/2) 



1 - 2a 



for a 7^ i. 



(21) 



(22) 



Furthermore, - give s m 2k. Figure [T] shows 
the plots of 7/_ as a function of k for r = 10 -4 and 
a= i ill 

u 2 ' 3 ' 4 ' 5 ' 

The following are some remarkable features of NSSs 
of the J-function potential (JTHJ) that distinguishes them 
from their linear counterpart. Notice that they hold ir- 
respective of the form of the nonlincarity profile f(\ip\). 

1. The condition for the creation of a NSS is highly sen- 
sitive to the value of a. 

2. Depending on the sign of 7/-, there are specific spec- 
tral gaps for NSSs. For example, as shown in Figure [T] 
for 7/_ > and a = 1/3, no NSS arises for the k/ir 
values in the intervals [1.5,2], [3,4], [4.5,5], [6,7], [7.5,8], 
[9,10], etc. In addition, the interval [0,0.5] is forbidden 
for all values of a whenever 7/_ > 0. For 7/_ < 0, NSSs 
reside on the spectral gaps of the case 7/- > 0. In par- 
ticular, LSSs are continuously related to the NSSs of the 
case 7/_ < 0. 

3. There is always a minimum value of I7/- 1 below which 
no NSS arises for k > 

for a = i. 

"nonlinearity threshold (NT)" . 

Suppose that the above model provides a description of 
an optical system consisting of a very thin planar slab of 



> ? if a ^ i and k > tt 
We will refer to this value of I7/-I as the 



2(l-2a) 



4 



k//r 



0.01 



-0.01 



0.01 



-0.01 



0.01 



-0.01 



0.01 



-0.01 











u 






' v. 






' — > 






n 

i 




n 
i i 




J 






I j 
u 






\ i 




r 






A 

D 


n 
i 




r 
i 






i 

j 


i 

u 






I 


1 
i 




r 
I 






\l I 






f 1 / 






J 


ji 


J 




V 


in 




< > 








n 


r 





k//r 

FIG. 1: (Color online) The location of the NSSs in the £-7/- 
plane for r = 10 -4 and a — ~, ~, |, | respectively from top 
to bottom. The vertical lines are asymptotes corresponding 
to the k/ir values violating the condition (|22[) . 



high-gain material. Because s is proportional to the gain 
coefficient @, we can adjust the value of s by controlling 
the pumping intensity. If \"/f~\ < NT, the system does 
not lase and the incident wave does not undergo a sub- 
stantial amplification regardless of how large s is. Now, 
suppose that \jf- \ > NT. Then as we increase s starting 
from zero, we find no amplification of the incident wave 
unless s reaches 2k\, where k\ is the smallest value of k 
such that (k, 7/-) corresponds to a NSS. Because we can 



use (fl~9|) to relate the values of 7/- and k, we can deter- 
mine one in terms of the other. For a Kerr nonlinearity, 
where 7/- = j\N^\ 2 , we can, in principle, employ this 
scheme to determine the frequency of the (incident) wave 
in terms of the amplitude of the transmitted wave. 

The rich structure depicted in Figure [1] suggests 
other potential applications of NSSs. For example, the 
parameter a, that signifies the center of the ^-function 
potential, can also be used as a control parameter in an 
experimental study of the above-mentioned frequency 
measurement scheme. Another possibility is to use 
independent frequency and intensity measurements 
together with the information about the location of 
NSSs to determine the coefficient of the Kerr and higher 
order nonlinearities of the medium. 

Concluding Remarks — In this letter we introduced the 
concept of a NSS for arbitrary confined nonlinearities and 
explored their properties for potentials having a compact 
support. In particular we explored in some detail NSSs 
of a complex (5-function potential and showed that they 
had a much richer structure than their linear counter- 
parts. Our results for this very simple model suggest, 
among other possibilities, a method for determining the 
frequency of an incident wave by performing an ampli- 
tude measurement. 

The results we report here may be viewed as a first 
step toward the study of the applications of nonlinear 
spectral singularities in various areas of physics. This 
might for example lead to the discovery of the analogs 
of threshold lasing and antilasing for nonlinear fields 
such as those encountered in acoustics, Bose-Einstein 
condensates, fluid mechanics, and even gravity. 
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